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First,  new  stability  tests  were  developed  for  multidimensional  recur- 
sive digital  filters  and  any  double-ended  n-dimensional  noncausal  linear 
processor  which  is  said  to  be  stable  if  its  impulse  response  decreases  ex- 
ponentially in  all  2-n  directions.  It  was  then  shown  that  the  impulse  response 
operator  for  a 2-D  discrete  Hilbert  transformer,  though  not  by  itself  sum- 
separable,  becomes  so  after  appropriate  classification.  Subsequently  it  was 
proved  that  the  multiplicative  complexity  of  computation  of  a 2-D  DHT  is  not-»^ 
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I 20.  Abstract  continued 

greater  than  twice  the  sura  of  multiplicative  complexities  of  two  1-D  DHT's. 
Subsequently,  the  1-D  matrix  Fade  approximation  problem  via  a three-term 
recursive  computation  scheme  was  tackled  as  a prelude  to  the  solution  of 
2-D  and  n-D  cases.  Specifically,  given  a l-D^^trix  power  series,  it  was 
shown.' that  a recurrence  relation  relates  the  , [^/M] , [^^>1^-])] 

order  Fade  approximants,  which  are  guaranteed  to^  exist  provided  a certain 
rank  condition  is  satisfied  by  characterizing  matrices  possessing  block- 
Hankel  structure.  Attention  to  stability,  algebraic  computational  com- 
plexity and  approximation  were  necessary  because  efficient  implementation 
of  stable  recursion  is  desired. 
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b.  STATEMENT  OF  WORK 

The  key  objective  of  the  research  has  been  to  make  fundamental  contribu- 
tions towards  the  design  of  multidimensional  recursive  digital  filters  and 
linear  processors.  The  major  emphasis  was  planned  for: 

(a)  development  of  new  multidimensional  stability  test  procedures  as 
well  as  stability  tests 

(b)  the  two-dimensional  discrete  Hilbert  transform  (DHT)  and  algebraic 
computational  complexity  aspects  for  Its  efficient  Implementation 
In  problems  of  2-D  stabilization  as  well  as  2-D  filtering,  and 

(c)  the  approximation  of  a prescribed  sequences  of  Impulse  responses 
that  could  characterize  either  a 2-D  single-input  single-output 
filter  or  a 1-D  multi- input  multi-output  (MIMO)  filter  as  a pre- 
lude to  the  approximation  of  multivariable  multidimensional  Impulse 
response  array  sequences. 

It  might  be  noticed  that  the  above  three  stated  objectives  are  closely 
linked  and  could  demand  attention  at  various  stages  of  a successful  design. 
Most  of  the  stability  procedures  developed  to-date  are  applicable  to  the 
bounded- input  bounded-output  (BIBO)  type  and  tend  to  become  computationally 
tedious,  especially  with  Increase  In  the  number  of  dimensions  of  the  filter. 
The  need  for  new  procedures,  therefore.  Is  not  over.  Also,  the  BIBO 
stability  criterion  may  not  be  sufficient  in  many  applications.  For  example, 
an  n-dlmenslonal  linear  processor  (noncausal,  in  general)  is  stable  if  its 
response  decreases  exponentially  in  all  2-n  directions.  This  can  be  trans- 
lated into  the  absence  of  zeros  of  the  denominator  polynomial  on  the  poly- 
dlsc  distinguished  botindary.  Thus,  the  potentially  trouble-causing  non- 
essential  singularities  of  the  second  kind  on  the  polydlsc  distinguished 
boundary  In  BIBO  stability  are  automatically  eliminated.  Applications  of 
this  type  of  stability  condition  Include  2-D  image  processing,  3-D  Image 
reconstruction  from  2-D  projections  and  real  time  Image  processing  of 
moving  objects.  Therefore,  a goal  of  the  research  has  been  to  develop  tests 
for  determination  of  the  presence  or  absence  of  multivariate  polynomial 
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zeros  on  the  unit  polydlsc  distinguished  boundary. 

After  application  of  the  stability  tests,  if  a filter  is  found  unstable, 
one  might  like  to  stabilize  the  unstable  filter  without  any  appreciable 
distortion  in  the  filter's  desired  amplitude  response  spectrum.  The  in- 
feasiblllty,  in  general  of  stabilizing  a 2-D  filter  using  extensions  of 
known  1-D  procedures  so  that  the  amplitude  response  exactly  matches  that 
of  an  unstable  filter  has  been  demonstrated.  In  the  absence  of,  as  yet, 
a general  2-D  stabilization  procedure,  it  appears  that  stabilization  via 
2-D  DBT,  whenever  applicable,  cannot  be  overlooked.  This  and  other 
potential  applications  including  2-D  lowpass  digital  filtering  provided 
the  motivation  for  planning  research  on  the  2-D  DHT.  The  goal  has  been 
to  study  the  special  properties,  if  any,  of  the  2-D  DHT  impulse  response 
operator  and  exploit  to  advantage  the  algebraic  structure  of  the  system 
in  order  to  be  able  to  implement  efficiently  the  2-D  discrete  Hilbert 
transformer,  kecent  results  on  the  minimal  multiplicative  complexity 
have  been  planned  for  use  in  the  present  application. 

The  approximation  of  functions  of  several  variables  is  probably  the 
branch  of  the  subject  offering  the  toughest  challenge.  Research  Into 
stability  and  stabilization  problems  were  planned  because  rational  approxl- 
mants  must  satisfy  stability  constraints  prior  to  recursive  Implementation. 
Though  much  of  the  one  variable  theory  can  be  extended  without  very  little 
modification,  fundamental  difficulties  often  occur  in  the  multidimensional 
arena.  For  example,  there  are  no  Chebyshev  sets  In  many  dimensions  so  that 
mlnlmax  approximation  is  not  necessarily  unique.  Non-uniqueness  does  not 
necessarily  disturb  one  in  practical  applications.  The  problem  of  unique- 
ness, however,  is  closely  connected  with  that  of  characterization,  and  the 
latter  is  of  much  greater  practical  importance.  Therefore,  multidimensional 
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■Inlnax  approxlmaclon  has  nostly  been  avoided  In  practice.  After  careful 
study  of  the  only  available  results  In  the  area  along  with  the  extension 
of  1-D  differential  correction  scheme  and  some  hybrid  approximation  schemes. 
It  was  decided  that  in  the  context  of  multidimensional  digital  filter  de- 
sign, development  of  recursive  schemes  for  lnq>lementlng  the  Fade  approxima- 
tion theory  might  be  the  most  appropriate.  Therefore,  research  Into  the 
theory  of  existence,  non-uniqueness  and  recursive  computation  of  2-D  scalar 
and  1-D  matrix  Fade  approxlmflnts  was  planned  as  a prelude  to  the  development 
of  the  corresponding  results  required  In  the  design  of  multivariable  multi- 
dimensional digital  filters  with  stability  constraints. 
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C.  STATUS  OF  RESEARCH  EFFORT 

The  overall  objective  of  the  research  has  been  to  contribute  towards 
the  approximation  problem  In  the  design  of  stable  multidimensional  recur- 
sive filters  and  linear  processors.  The  key  factors  considered  In  this 
research  towards  the  fulfillment  of  the  final  objectives  are  development 
of  new  multidimensional  stability  tests,  the  approximation  of  a power 
series  in  several  variables  by  rational  approxlmants  and  the  algebraic 
computational  complexity  problems  that  must  be  attended  to  for  efficient 
Implementation . 

1.  Multidimensional  Stability  Tests:  Over  the  last  few  years  multidimensional 

systems  research  has  progressed  along  different  directions  as  described  In 

[1]>  [2].  The  topic  of  stability  tests  for  multidimensional  filters  has 

been  the  subject  of  Intensive  research  over  the  past  several  years  and  an 

up-to-date  exposition  of  research  on  the  subject  Is  available  in  [3].  Most  of 

the  tests  have  centered  around  bounded- input  bounded-output  type  of  stability, 

while  an  n-dimenslonal  linear  processor  Is  stable  If  Its  impulse  response 

decreases  exponentially  In  all  2-n  directions.  It  can  be  shown  that  every 

rational  processor  transfer  function  has  a unique  stable  expansion  If  the 

denominator  polynomial  B(z. , z-,  ...,z  ) in  the  complex  variables  z, , z. , . . .,  z 

J-  z n 1 2 1 

satisfies  the  condition  given  in  (1). 

*2 *n)  ^ 0.  on  Iz^l  - 1*2!“  • • • “ l*nl  " ^ <1) 

Applications  include  2-D  Image  processing,  3-D  image  reconstruction  from  2-D 
projections  and  cross-sections,  real  time  image  processing  of  moving  objects, 
repetitive  Infinite  networks  and  possibly  in  tests  for  existence  of  Green's 
function  with  vanishing  infinite  boundary  conditions  in  linear  partial 
difference  equations.  The  present  research  contributed  to  the  solution  for 
the  teat  of  condition  in  (1),  which  did  not  exist  except  for  n ■ 1 case. 
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IWo  Indirect  and  one  direct  procedure  to  test  a prescribed  multivariate 
polynomial  for  zeros  on  a polydlsc  distinguished  boundary  substantiate 
constructively  the  existence  results  from  elementary  decision  algebra  [4]. 

The  exact  methods  lmplementable,lf  necessary,  with  arbitrary  precision, 
complement  the  numerical  method  Involving  DFT  computations.  In  many  cases, 
the  latter  method  can  give  satisfactory  result  with  speed  while  In  some 
cases  Involving  high  coefficient  sensitivities  (not  unconmion  in  multi- 
dimensional problems)  the  procedures  given  In  the  paper  [5],  resulting 
from  this  research,  are  more  suited.  In  this  paper,  the  first  Indirect  test 
formulation  relates  the  present  test  to  an  already  developed  global  positivity 
test  on  a polynomial  in  several  real  variables  [6]  while  the  second  relates 
the  present  test  to  an  already  developed  local  positivity  test  package  [7], 

[8].  A direct  test  Is  then  formulated  for  the  problem  under  consideration. 
Detailed  discussions  are  given  for  the  n >•  2 case  with  adequate  attention 
to  the  various  singular  cases  chat  might  occur  In  the  implementation  of 
the  procedure.  The  feasibility  of  implementing  the  test  via  rational  opera- 
tions is  shown  for  the  n - 2 case.  Parallels  for  the  n > 2 case  are  briefly 
given,  and  the  details  with  Illustrative  example  for  the  direct  test  along 
with  Its  limitations  when  n > 2,  Is  available  In  [9].  The  Important  research 
results  obtained  In  the  area  have  been  reported  in: 

N.  K.  Bose  and  S.  Basu,  "Test  for  polynomial  zeros  on  polydlsc  distin- 
guished boundary,"  IEEE  Transactions  on  Circuits  and  Systems  Special  Issue 
on  '^Mathematical  Foundations  of  System  Theory,"  Vol.,25,  Sept.,  1978,  pp. 684-693. 

Another  n-dimenslonal  filter  bounded- Input  bounded-output  stability  test 
was  developed  during  the  course  of  the  research  as  a follow-up  of  a new 
stability  test  and  Its  Implementation  for  two-dimensional  filters,  advanced 
by  the  principal  investigator  earlier  [10].  It  was  shown  how  after  a single 
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dimensional  stability  test,  the  final  result  obtained  for  the  n-D  case  re- 
quires the  testing  of  a single  variable,  two  variable,  ...  (n-1)  variable 
polynomials  (in  real  variables)  for  local  positivity  In  -1  ^ ^ 1, 

-1  ^ ^ 1,  1 ■ 1,2,  ...  * “1  Xji^  ^ 1,  1 ■ 1,2,  ...  , n-1.  These  poly- 

nomials can  be  conveniently  generated  in  the  sense  that  all  remaining  poly- 
nomials can  be  obtained  by  specializing  to  zero  relevant  variables  at 
individual  stages  In  the  recursive  scheme  used  In  the  computation  of  the  (n-1) 
variable  polynomial.  The  research  results  obtained  in  this  connection  were 
published  In: 

M.  K.  Bose,  "Implementation  of  a new  stability  test  for  n-D  filters," 

IEEE  Transactions  on  Acoustics,  Speech  and  Signal  Processing,  Vol. , 27, 
Febrxiary,  1979,  pp.  1-4. 

2.  2-D  DHT  and  Computational  Complexity  Aspects;  It  was  shown  that  the 

2-D  discrete  Hilbert  transform  (DHT)  operator  may  be  expressible  as: 

H(i^,  ij)  - jj-J—  [{1  - (-1)  + {1  - (-1)^}I2  - - (-1) 
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where 


ill 

N, 


ill 


j - , Wj^  - e ^ , w 


In 


*1  , ,,h  *2  , ,,^2 

• "l  - . I;  • “2  - 


V * 1 
1 


In  the  standard  vector  form  representation, 
vec  (PI)  - H vec  {PR}  (or  Y - H X) 
of  the  2-D  convolution. 


Ni-1 


N2-I 


Pldl’  " k^-o  k^-O  PR(kj^,  k^)  H(ij^  - k^^,  1^  - k^) 

for  Ij^  - 0,1,2,  . . . , - 1,  1^  - 0,1,2 - 1, 


t 
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the  IL  ® block-circulant  structure.  Each  of  the 

N2  ^ element  matrices  In  ^ possesses  a clrculant-llke  structure.  This 
algrebralc  structure  was  exploited  to  arrive  at  the  main  results  summarized 
In  lemma  '<  and  theorem  1 below. 

Lemma  1.  The  2-D  DHT  of  an  Input  array  can  be  computed  using 

only, 

oCNj^,  N^)  - 2o(Nj^)  + 2a(N2) 

multiplications,  where  a(N)  Is  the  number  of  multiplications  required  In 
the  computation  of  the  1-D  DHT  of  a N-polnt  Input  vector. 

Theorem  1.  The  2— D DHT  of  an  Input  array  can  be  computed  using 

only 

a(N^,  N^)  - A(N^  + N^)  - 2(4(N^)  + ♦(N^)  ) 
multiplications  (provided  multiplication  by  any  fixed  element  of  the  field 
Is  not  counted  as  1 multiplication),  where  i^(N)  Is  the  number  of  distinct 
Irreducible  factors  over  G of  the  polynomial  in  u,  namely  u^  - 1. 

In  the  absence  of,  as  yet,  a general  2-D  stabilization  procedure,  it 
appears  that  stabilization  via  2-D  DHT,  whenever  applicable,  cannot  be 
overlooked.  This  and  other  potential  applications  Including  2-D  lowpass 
digital  filtering  provided  motivation  towards  the  obtaining  of  results  on 
the  two-dimensional  DHT  and  computational  complexity  aspects  In  Its  Im- 
plementation as  an  Integral  part  of  the  topic  of  the  present  research. 

The  details  of  the  research  Just  outlined  will  appear  In: 

N.  K.  Bose  and  K.  A.  Prabhu,  "2-D  DHT  and  computational  complexity 
aspects  In  Its  Implementation,"  accepted  for  publication  In  IEEE  Trans- 
actions Acoustics,  Speech  and  Signal  Processing,  August,  1979. 

In  view  of  the  Important  role  of  finite  field  theory  In  the  efficient 
implementation  of  developed  algorithms,  attention  was  directed  towards  some 
other  aspects,  which  show  strong  potentialities  for  future  engineering 
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applications  and  offer  scope  for  further  investigations  and  developments. 

In  texts,  a procedure  to  obtain  the  number  of  distinct  Irreducible  monlc 
polynomials  of  prescribed  degree  m in  a single  variable  x^^,  having  coeffi- 
cients over  a finite  field  F of  q elements,  has  been  given.  This  he.>  been 
done  through  the  use  of  the  concepts  of  generating  functions  and  Moeblus 
inversion.  The  result  has  been  of  use  In  multiple  error-correcting  code 
construction  as  well  as  in  other  problems  of  algebraic  coding  theory.  In 
the  present  research  a procedure  to  determine  the  number  of  distinct  ir- 
reducible polynomials  of  prescribed  degrees,  m, , m , . . . , m in  the 

1 z n 

Independents  variables  x , x , . . . , x , respectively,  having  coefficients 
over  a finite  field  of  q elements  was  obtained.  It  was  seen  that  certain 
interesting  distinction  from  known  results  occur  in  the  resolution  of  this 
problem — a fact  often  encountered  in  the  transition  from  the  single 
dimensional  to  the  multidimensional  case.  Though  the  results  arrived  at 
have  their  validity  for  an  arbitrary  n,  in  the  research  conducted,  for  the  sake 
of  brevity  in  exposition  attention  has  been  especially  concentrated  on  the 
n ■ 2 or  the  2-D  case  as  reported  in  the  preprint: 

K.  A.  Prabhu  and  N.  K.  Bose,  "Number  of  irreducible  q-ary  polynomials 
in  several  variables  with  prescribed  degrees,"  accepted  for  publication  in 
IEEE  Transactions  of  Circuits  and  Systems. 

3.  Approximation  theory  results:  The  1-D  scalar  Fade  approximation  theory 
has  been  widely  used  in  problem  of  theoretical  physics  [11],  [12],  numerical 
analysis  [13],  and  electrical  engineering  especially  because  of  its  close 
link  to  the  Prony  method  [14]  and  the  theory  of  continued  fractions  [15]. 

The  necessity  for  generalizing  the  1-D  Fade  scalar  theory  to  the  1-D(m  x n) 
matrix  case  was  discussed  in  the  setting  of  the  minimal  partial  realization 
problem  [16].  The  feasibility  for  recursive  computation  of  the  1-D  matrix 
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Fade  approxlmants  in  the  context  of  minimal  partial  realization  theory  was 


shown  in  [17].  The  matrix  Fade  approximation  theory  In  the  setting  of  a 


noncommutatlve  algebra  with  a view  towards  application  in  scattering  physics. 


Is  also  available  in  [1,  ch.  21],  [12,  pp.  19-44],  Hore  applications  of 


the  1-D  matrix  Fade  theory  can  occur  In  multiport  network  synthesis  and  In 


the  design  of  multi- Input  multi-output  digital  filters.  The  present 


research  provided  an  efficient  recursive  computational  scheme  for  the 


1-D  matrix  Fade  approxlmants  as  a prelude  to  2-D  and  n-D  results  required 


In,  among  other  problems,  in  the  design  of  multivariable  multidimensional 


digital  filters.  The  main  contribution  of  the  present  research  was  the 


demonstration  of  the  existence  of  a recurrence  relation,  using  which  the 


[ /M+1]  order  1-D  matrix  Fade  approximant  can  be  computed  from  the 


[^/M]  and  [^  ^/M-1]  order  1-D  matrix  Fade  approxlmants  to  a prescribed 


matrix  power  series  T(z)  of  order  m x n.  It  Is  also  possible  to  obtain  the 


matrix  polynomial  B , (z)  associated  with  the  [ /Mfl]  order  approximant 
MtI 


by  solving  a system  of  linear  equations,  characterized  by  a block  Hankel 


matrix  of  order  {(M+l)m}  x {(M+l)n}.  In  [18],  for  the  m • n ■ 1 case  a 


solution  scheme  requiring  0(M  ) operations  has  been  advanced  and  the 


feasibility  for  extension  to  the  more  general  block  matrix  case  has  only 


been  cited.  Therefore,  computational  complexity  comparison  between  the  two 


schemes  can  be  carried  out  as  a separate  project.  The  advantage  of  the 


technique  presented  here  appears  to  lie  In  the  use  of  the  recurrence 


relation  established,  where  the  order  of  matrices,  whose  Inverses  or  right 


Inverses  have  to  be  computed.  Is  Independent  of  the  order  of  approxlmants. 


The  details  of  these  results  are  available  in: 


N.  K.  Bose  and  S.  Basu,  "Theory  and  recursive  computation  of  1-D  matrix 


Fade  approxlmants,"  Froceedlngs  of  1978  Conference  on  Dec.  and  Control, 


pp.  653-657.  Also  an  article  has  been  completed  for  submission  to  a Journal. 
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tional complexity  aspects  in  Its  Implementation,"  IEEE  Transactions  of 
Acoustics,  Speech  and  Signal  Processing,  accepted  for  publication  in 
August,  1979.  • 

K.  A.  Prabhu  and  N.  K.  Bose,  "Number  of  Irreducible  q-ary  polynomials 
In  several  variables  with  prescribed  degrees,"  IEEE  Transactions  of 
Circuits  and  Systems,  accepted  for  publication  in  1979. 
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N.  K.  Bose,  ed. , "Multidimensional  Systems;  Theory  and  Applications." 
IEEE  Press,  N.Y.,  1979. 


e.  RESEARCH  PERSONNEL 


Dr.  N.  K.  Bose,  Professor  of  Electrical  Engineering  and  Mathematics 
served  as  the  Principal  Investigator  throughout  the  period  of  report 
extending  from  January  1,  1978,  to  January  31,  1979.  Worked  and 
guided  work  on  problems  of  multidimensional  stability,  stabilization, 
approximation  and  algebraic  computational  complexity  in  implementation. 
Mr.  S.  Basu,  a Graduate  Student  Researcher,  worked  on  research  supported 
by  the  grant  for  the  complete  period.  He  completed  a M.S.  thesis 
entitled,  "Theory  and  applications  of  a direct  test  procedure  for  poly- 
nomial zeros  on  polydisc  distinguished  boundary."  at  the  University  of 
Pittsburgh  In  April,  1978.  Subsequent  to  that,  he  worked  on  the 
approximation  problem  in  the  design  of  multidimensional  digital  filters. 
Mr.  K.  A.  Prabhu,  worked  for  five  months  as  a Research  Fellow  and  one 
month  as  a Graduate  Student  Researcher.  His  efforts  were  concentrated 
primarily  towards  the  algebraical  computational  complexity  aspects  In 
the  implementation  of  stability  and  approximation  algorithms  occurring 
in  the  design  of  multidimensional  filters.  Approximately  two  man- 
months  of  support  were  also  provided  by  Graduate  Student  Researchers, 

Mr.  Wen  and  Mr.  John  W.  Loney.  Mr.  Wen  began  to  look  Into  the  efficient 
implementation  of  2-D  algorithms  for  signal  processing  while  Mr.  John  U. 
Loney  participated  for  a month  on  discussions  involving  scopes  for 
multidimensional  minlmax  approximation  which  was,  however,  discarded 
in  favor  of  Pade  approximants  with  stability  constraints  as  more 
suited  for  the  type  of  research  undertaken  here. 
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f.  INTERACTIONS  OF  PRINCIPAL  INVESTIGATOR.  DR.  N.  K,  BOSE 

1.  Invited  to  be  the  Organizer  and  Session  Chairman  of  a Special  Session 
on  "Multidimensional  Systems,"  at  the  1978  International  Symposium  on 
Circuits  and  Systems  held  at  New  York  City,  May  16-19,  1978.  Also 
presented  a paper  entitled: 

"Test  for  polynomial  zeros  on  polydlsc  distinguished  boundary  and 
consequences," 

at  the  Symposium. 

2.  Invited  for  collaborative  research  on  multidimensional  wave  digital 
filters  for  one  month  at  Ruhr  Unlversltat,  West  Germany  In  1978. 

Visit  was  sponsored  by  Deutscher  Akademlscher  Austauschdlenst.  Took 
part  In  numerous  technical  discussions  with  Professor  Fettwels  and  his 
assistants  - primarily  Herr  Klaus  Meerkotter  and  also,  to  some  extent, 
Herr  Lennarz  and  Herr  Llnnenberg.  The  topics  of  discussion  centered 
primarily  around  various  problems  associated  with  the  synthesis  of 
two-dimensional  wave  digital  filters,  as  well  as  the  multidimensional 
stability  tests  recently  developed  by  the  author,  and  the  importance 

of  approximation  In  the  context  of  2-D  filter  design. 

3.  Has  an  organizer  and  participant  In  the  National  Science  Foundation 
sponsored  Regional  Conference  on  Algebraic  Complexity  of  Computation 
and  Its  Applications  to  Problems  In  Engineering  and  Computer  Sciences, 
held  at  the  University  of  Pittsburgh  during  Aug.  21-25,  1978.  Intro- 
duced here  the  properties  of  2-D  Hilbert  transformer  and  subsequently 
employed  aspects  of  computational  complexity  theory,  presented  by 

S.  Wlnograd,  to  efficiently  Implement  the  2-D  discrete  Hilbert  trans- 
form In  2-D  stabilization  and  2-D  digital  filtering  problems. 
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Invited  to  present  a paper  In  the  Special  Session  on  "Hultldlmenslonal 
Systems"  organized  by  Dr.  M.  Korf  of  Stanford  Uhlverslty  and  Dr.  B.  Levy 
I of  the  Massachuaetts  Institute  of  Technology  at  the  Conference  on 

^ Decision  and  Control  held  at  San  Diego,  Jan.,  10-12,  1979.  The 

I Invited  talk  presented  at  the  Conference  was  entitled,  "Multidimensional 

Systems  Theory:  Present  state  as  an  indicator  of  future  prospects," 
end  a portion  of  the  talk  was  published  In  the  Conference  Proceedings 
os,  "Multidimensional  Systems  Theory:  1— D Matrix  Pade  Approxlmants." 

5.  The  paper  "Sum  and  product  separabilities  of  multivariable  functions 
and  applications,"  by  S.  Chakrabarti,  N.  K.  Bose  and  S.  K.  Mltra  was 
Included  in  the  book  entitled,  Two-Dimensional  Digital  Signal  Proces- 
elng.  published  by  Academic  Press  in  1978  as  volume  20  in  Benchmark 
Papers  in  Electrical  Engineering  and  Computer  Science. 
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g.  NEW  RESULTS  AND  SPECIFIC  APPLICATIONS  RESULTING 
FROM  RESEARCH  EFFORT 


1.  The  tvo  Indirect  and  one  direct  test  procedures  for  seros  on  the  dis- 
tinguished boundary  of  the  polydlsc  provide  tests  for  stability  of  a 
n-dlmenslonal  linear  processor  which  Is  required  to  have  Its  Impulse 
response  decrease  exponentially  In  all  2-n  directions.  Applications 
Include  2-D  image  processing,  3-D  image  reconstruction  from  2-D  pro- 
jections, real  time  image  processing  of  objects  or  targets  In  motion, 
repetitive  Infinite  networks  and  tests  for  existence  of  Green's  function 
with  vanishing  infinite  boundary  conditions  In  linear  partial  difference 
equations.  The  Indirect  tests  dwell  on  the  multivariate  polynomial 
global  positivity  and  local  positivity  test  procedures  previously 
developed  by  the  principal  investigator. 

2.  The  new  BIBO  stability  test  procedure  for  multidimensional  recursive 
filters  has  a natural  potential  for  application.  The  efficiency  in 
Implementation  and,  consequently,  the  popularity  of  the  procedure  will 
increase  with  the  improvement  of  multivariate  polynomial  local  positivity 
test  methods. 

3.  The  demonstration  of  the  fact  that  the  multiplicative  complexity  for 
implementation  of  a 2-D  DHT  is  not  greater  than  twice  the  multiplicative 
complexities  of  only  two  1-D  DHT's  will  be  a strong  motivation  behind 
the  use  of  2-D  DHT  In  stabilization  and  low-pass  digital  filtering 
problems. 

4.  The  development  of  a recursive  enumeration  procedure  for  the  number 
of  Irreducible  q-ary  polynomials  In  several  variables  with  prescribed 
degrees  is  a fundamental  contribution  In  multidimensional  signal 
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^ procesalng  over  finite  fields.  Scopes  for  application  are  Included 

I 

in  signal  processing,  system  Identification,  pseudo-random  noise 
generation,  q-ary  array  synthesis  and  multidimensional  codeword  gen- 

I 

I aratxon. 

5.  Finally,  the  development  of  a three-term  recursion  procedure  for 

I 

I 

computation  of  1-D  matrix  Fade  approxlmants  Is  likely  to  be  of  use 
In  several  areas  besides  MIMO  digital  filtering.  This  fundamental 
contribution  can  be  directly  applied  to  the  minimal  partial  realiza- 
tion problem.  Identification  of  linear  dynamical  systems  from  time- 
domain  measurements,  the  scattering  problem  in  high  energy  theoretical 
ptqfslcs,  and  possibly  in  the  synthesis  of  multiport  distributed  net- 
works by  lumped  approxlmants. 


